Non-Gaussianity of a single scalar field in general covariant Hofava-Lifshitz gravity 
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In this paper, we study non-Gaussianity generated by a single scalar field in slow-roll inflation in 
the framework of the non-relativistic general covariant Hofava-Lifshitz theory of gravity with the 
project ability condition and an arbitrary coupling constant A, where A characterizes the deviation 
of the theory from general relativity (GR) in the infrared. We find that the leading effect of self- 
interaction, in contrary to the case of minimal scenario of GR, is in general of the order a n e 3 ^ 2 , where 
e is a slow-roll parameter, and a n (n = 3, 5) are the dimensionless coupling coefficients of the six- 
order operators of the Lifshitz scalar, and have no contributions to power spectra and indices of both 
scalar and tensor. The bispectrum, comparing with the standard one given in GR, is enhanced, and 
gives rise to a large value of the nonlinearity parameter /nl • We study how the modified dispersion 
relation with high order moment terms affects the evaluation of the mode function and in turn 
the bispectrum, and show explicitly that the mode function takes various asymptotic forms during 
different periods of its evolution. In particular, we find that it is in general of superpositions of 
oscillatory functions, instead of plane waves like in the minimal scenario of GR. This results in a 
large enhancement of the folded shape in the bispectrum. 



I. INTRODUCTION 

The Horava-Lifshitz (HL) theory of quantum gravity, 
proposed recently by Horava [1], motivated by the Lif- 
shitz scalar field theory in solid state physics [2], has 
attracted a great deal of attention, due to its several re- 
markable features in cosmology as well as some challeng- 
ing questions, such as instability, ghost, strong coupling, 
and different speeds Q. To resolve these questions, vari- 
ous models have been proposed, along two fundamental], 
different lines, one with the projectability condition 

N = N(t), 



it ally 

m, 
(i.i) 



and the other without it 0, [TO], where N denotes the 
lapse function in the Arnowitt, Deser and Misner de- 
compositions [Til] - I n particular, Horava and Melby- 
Thompson (HMT) proposed to enlarge the foliation- 
preserving diffeomorphisms of the original model, 



St 



-/(*), ^ = -C(t,x), 



(1.2) 



often denoted by Diff(M, J 7 ), to include a local U(l) 
symmetry, so that the reformulated theory has the sym- 
metry 0, 



U(l) k Diff(M, T). 



(1.3) 



With such an enlarged symmetry, the spin-0 gravitons, 
which appear in the original model of the HL theory 
[lj], are eliminated 0, 0, and as a result, all the prob- 
lems related to them, including instability, ghost, and 
strong coupling, are resolved automatically. This was 
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initially done with A = 1, where A characterizes the de- 
viation of the theory from general relativity (GR) in the 
infrared, as one can see from Eqs. ([2.1|) and (|2.2p given 
below. Soon, it was generalized to the case with any 
A [7[, in which it was shown that the spin-0 gravitons 
are also eliminated 0, 0], so that the above mentioned 
problems are resolved in the gravitational sector. In the 
matter sector, the strong coupling problem, first noted 
in [8[, can be solved by introducing a mass M* so that 
M* < A w , where M* denotes the suppression energy of 
high order operators, and A w the would-be energy scale, 
above which matter becomes strongly coupled [12j, sim- 
ilar to the non-projectability case without the enlarged 
symmetry [13|- The consistence of this model with so- 
lar system tests was investigated recently Q ? an d found 
that it is consistent with all such tests, provided that the 
gauge field and Newtonian prepotential are part of the 
metric, as proposed originally by HMT [5[- 

In this paper, we shall work within the HMT frame- 
work of the HL theory with the projectability condition 
[5-0. For the current status of the models proposed in 
[4 , [qI, |To| , we refer readers to [l4[ . 

Since the HL theory differs from GR significantly in 
high energies, in this paper we study another important 
issue - the primordial non-Gaussianity in the HL theory. 
Such a problem was studied previously, one without the 
projectability condition [l5[, and the other with it pij 
but in the curvaton scenario [ItJ, and many interesting 
results were obtained. However, our study presented in 
this paper are different from these at least in two aspects: 
(i) we shall study the problem in the framework of the 
non-relativistic general covariant Horava-Lifshitz theory 
of gravity with the projectability condition [5-8], as men- 
tioned above, in which the degree of gauge freedom is the 
same as that in GR; and (ii) we shall investigate the prob- 
lem in the inflationary scenario for a single scalar field 



2 



with slow-roll conditions 1 . Because of these differences, 
our results are also significantly different from theirs. 

To study non-Gauusianity, various techniques have 
been developed. In particular, the in-in formalism, de- 
veloped initially by Schwinger some decades ago [20| , be- 
comes standard, after it was first explored by Madalcena 
in his pioneer work [2l[ in calculating the high-order cor- 
relators for cosmological perturbations, and then further 
developed by Weinberg in his seminal paper [22]. For 
reviews, see [23|, [24[. In this paper, we will follow this 
approach. 

Non-Gauusianities have attracted lot of attention re- 
cently and been studied extensively in various models 
[23l . Ejij . mainly because they could be well within the 
range of detection of the current Planck satellite [25| and 
the forthcoming experiments, such as the CMBPol mis- 
sion [26j. Among various models, the one with a canoni- 
cally coupled single scalar field in the framework of GR, 
though the simplest, gives elegant predictions and fits 
extremely well with current observations [27j. It pre- 
dicts that the quantum fluctuations, which are respon- 
sible for generating the primordial perturbations and in 
turn the anisotropics in CMB and inhomogeneity in the 
large scale structure, are largely Gaussian. The effect of 
the non-linearity in the primordial perturbations is an 
order smaller than that of the power spectrum and is 
beyond our current detect abilities. 

However, Holman and Tolley recently argued that, 
while a change in high energy physics during or before 
inflation gives mild modifications to the power spectrum 
and index, the non-Gaussianity, evaluated by higher or- 
der correlation functions of the perturbations, is much 
more sensitive to the new physics in the ultraviolet (UV) 
(28| . Therefore, it is very interesting to study the non- 
Gaussianity in the scalar primordial perturbations of the 
HL theory, where the dispersion relations are quite dif- 
ferent from the standard one. 

The rest of the paper is organized as follows. In Sec- 
tion II we give a brief review of the non-relativistic gen- 
eral covariant Horava-Lifshitz theory of gravity with the 
projectability condition and an arbitrary coupling con- 
stant A. The self-interaction Hamiltonian and the lead- 
ing order terms are analyzed in Section III. Section IV 
discusses the modified dispersion relation and its effect 
on the mode function and the shapes of bispectrum. We 
also compare our results with those presented in [l6| . 
Finally, in Section V we summarize our results and the 
assumptions made along the way of analysis. Three ap- 
pendices are also included, where in Appendix A the lin- 
ear perturbations and the corresponding field equations 
are presented, while in Appendix B, the cubic part of 
the total action is calculated, and given explicitly. In 



Appendix C, the matching across the boundaries among 
the three different regions of the time revolution of the 
mode function [cf. Fig. [T] are given. 



II. GENERAL COVARIANT HL GRAVITY 
WITH PROJECTABILITY CONDITION 



The action of the general covariant HL theory of grav- 
ity with the projectability condition can be written as 



S = J dt(PxNy/g(c K -Cv + C<p+£A + C 

dtd 3 xN^/g£ M , (2.1) 



where g = det(g^), G is the Newtonian constant, and 

C K = K id K ij -\K\ 

C v = C 2 g +giR+^ {g 2 R 2 + gsR^R lJ ) 



~ (g^R 3 + g 5 R RijR ij + geR^R?) 
~ [g 7 RAR + g 8 (V t R jk ) (V^)] , 



r = 

'-'in 



Ca = ^(2A g -R), 

C x = (1- \)\(A<p) 2 + 2KAip 



(2.2) 



Here A = g^ViVj, A^ is a coupling constant, the Ricci 
en 

metric g^ , and 



and Riemann tensors Rij and R l j kl all refer to the 3 



Kij = ^(-gij + ^iNj+VjNi), 
Gij Rij - -gijR + k g 9ij. 



(2.3) 



The coupling constants g s (s = 0,1,2, ...8) are all di- 
mensionless, and 



A = 



32ttG' 



(2.4) 



1 To solve the horizon problem and obtain almost scale-invariant 
perturbations in the HL theory, inflation is not necessary, as first 
noted in [l8|. But, to solve othe problems, such as the relics of 
topological defects, including monopoles, it is still needed [l9j . 



is the cosmological constant. The relativistic limit in the 
IR, on the other hand, requires, 



Ai = -i. 



(2.5) 



3 



Cm is the Lagrangian of matter fields, and for a scalar 
field x, it is given by [HJ S3, 



Cm = 4°> + 4 A '*>, 



4 0) 



V 



V, 



/(A) 
2N 2 

v (x) + ( $ + Vx ( X ) ) (v x ) 2 + v 2 ( X ) v'i 



(x - N'ViX 



1 



+v 3 (x) v\ + v 4 (x) r 2 + v 5 (x) (Vx) 2 p 2 

+V6(x)^i^2, (2.6) 
^4 - A r 

ci (x)Ax + c 2 (x) (Vx)' 



AT 

iV 

/ 

2 



(x-^V;x) (vV)(v fc x) 



(vV) (v fe x) 



(2.7) 



with C\ (x) , C2 (x) 5 ^(x) and V^(x) being arbitrary func- 
tions of x, and 



= A"X- 



(2.8) 



For detail, we refer readers to [8|, LL9j. In the following 
we shall use directly the symbols and conversions from 
these papers, without further explanations. 



III. NON-GAUSSIANITIES 

The homogeneous and isotropic flat universe is de- 
scribed by, 

N = a(rj), gij =a 2 (rj)6ij, Ni = = <p, i = i(ry), 

n (3 - 1} 

where as in [19| , we use symbols with bars to denote the 
quantities of the background in the (£, ^-coordinates, 
and the ones with hats to denote those in the (rj,x 1 )- 
coordinates, where the conformal time r] is defined as 
r] = J dt/a(t). The linear perturbations, given by 

SN a<j>, SNi a 2 B yi , 

Sgij = -2a 2 {^5ij - E yij ), 

A = A + aSA, ip = (p + Sep, 

X = X + ^X, (3.2) 

were studied in detail in [l9| , and shown explicitly that a 
master equation exists for a single scalar field, with the 
gauge choice, 



= = E = Sip. 



(3.3) 



For the sake of reader's convenience, we summarize the 
main equations obtained in [l9[ in Appendix A of this 
paper. 

To generalize the above linear perturbations to the 
nonlinear case, in this paper we consider the perturba- 
tions, given by 

gij = aV%, Ni=a 2 B^ N = N, 

A = A + a5A, <p = 0, x = X + ^X- (3.4) 



Clearly, to first order, they reduce to the ones given by 
Eq. (|3.2p if one identifies £ as £ = —ip. Substituting 
Eq. (j3.4|) into the total action (j2.1|) . we find that its cubic 
part is given by Eq. (jB.l|) in Appendix B. From the linear 
order constraint and field equations given by Eq. (|A.10|) - 
(jA.llj) in Appendix A, we can express the terms B and 
5 A in terms of £. Then, substituting them into the cubic 
action (jB.l|h we find that it can be cast in the schematic 
form, 

r 3 

S = J2 { d2k ^ m ' ^X (3 " m) }. (3.5) 

To find the leading order terms in the self-interaction, 
let us first note that [l9| . 



Sx = h1Z(x e 1/2 M pl ft, 



(3.6) 



where 7Z is the comoving curvature perturbations, de- 
fined explicitly in [19]. From Eq. (|A.12|h we also have 



2M p V 



(3.7) 



Assuming that c\ ~ M* <C M v \ [19|, we find that ( <C 
Sx <C TZ. This implies that to find the leading order 
terms, it suffices to look for terms which are of cubic order 
of Sx [ m — i n (|3.5p ]. since all terms in lower orders of 
Sx are of higher orders of £, hence further suppressed by 
factors of M*/M v \. With this as our guideline, it can 
be shown that only six terms are left for considerations. 
One from the part *5 X |^ in Eq. (|B.3|) , identified as 



^"'(X^X 3 . 



(3.8) 



However, since it is the third-order derivative term of the 
potential, one can immediately ignore it, as the slow-roll 
conditions require \V"'\ <C 1. The other five are from 
S x \^ and are identified as, 



a z a 4 a 4 



(d 2 5 X y,^(d 4 sx) (d k 6 X y. 



(3-9) 



Out of the five, three are proportional to derivatives of 
the coupling functions V± , V2 , and V4 . They all appear in 
the linear perturbations, and it was assumed that their 
derivatives with respect to x vanishes [l9j |. To be con- 
sistent with it, in this paper we keep this assumption. 
Hence, we are finally left only with two terms, that are 
proportional to V3 and V5, 



dr]d 3 x 



h 3 { V 3 (d 2 K)\v 5 (d 4 K) (d k Kf \, (3.10) 



where h is defined in Eq. (jA.16|) . In contrast to the 
minimum scenario in GR, which predicts that the self- 
interaction should be of the order of e 2 , in the current 



case, the leading order is of a n h 3 ~ d n e 3 / 2 , where a n are 
dimensionless parameters, defined by V n = a n /M^ (n = 
3,5). 

Despite the fact that these two terms are similar to 
the OL2 and as terms given in [16|, a key difference, how- 
ever, exists. In [lq , the authors worked in the framework 
of curvatons [17| , and inflation was not necessary to pro- 
duce the scale- invariant power spectrum. As a result, the 
time of interest was assumed to be the period in which we 
have H ^> M*. Thus, the quantization of their Lifshitz 
scalar (/>, which is responsible for generating the primor- 
dial curvature perturbations, can be carried out as [161 ] . 



0(x,t) 



5 

A* 




. / d 3 ke~ ikx 


(2tt) 3 J 






" . fc 3 / 


fc3/2 eX P 


. % M*J 



(it 
^2 



(3.11) 



The mode freezes after it leaves the sound horizon 
(i^M 3 ) -1 / 4 , which is much smaller than the Hubble hori- 
zon in the inflation scenario, and gives rise to a power 
spectrum, 



and find that to the leading order in H\ we have 

{n kl {t)n^{t)^At)) 

dt' ( [#| (3 ) (f), n kl (t)K k2 (t)n k3 (t) 

^3 r3 



ih 6 (27T)- 5 S 6 (ki + k 2 + k 3 ) 



6V 3 



^1^2 ^3 



v 5 



(^1 + ^2 "I" ^f) 



- {k\k\ + k\k\ + cyclic) 
x - c.c.j, 



(3.15) 



where 



tf(f;f) = r^(t)^ 2 (t)4(t) 
rt dt' 

^r fa (fK(fK(0- (3-16) 

r&(£) is the mode function for 7£ 2 . In writing down 
(j33SJ), we assumed that d(V 3 ^)/dt = V^x = 0. We'll 
see that, once this assumption is relaxed, it will generate 
more interesting features in the shapes of the bispectrum. 



M 2 



(2tt)V 



(3.12) 



This is quite different from the expression 
pGR _ hA 



(2^) 2 ^ 2 ' 



(3.13) 



obtained with inflation in GR. In contrast, slow-roll in- 
flation is required in our current model and the power 
spectrum was found to resemble the GR one [l9|, given 
by 



(3.14) 



where e HL H 2 /M% <C 1. 

With the self-interaction at hand, and the relation 
$1(3) = — / dtH\(s\, we perform the calculations of the 
bispectrum by using the in-in formalism [2l|, [22| , 



Qit) 



T exp(z jf Hi{t')dt'^j 



Q(t) 



T exp 



(-if Hi(t')dt' 

V Jto 



IV. MODIFIED DISPERSION RELATIONS AND 
THE SHAPES OF THE BISPECTRUM 

To study the shapes of the bispectrum, in this section 
we first consider the time evolution of the mode function 
in the (quasi-)de Sitter background. 



A. Time evolution of the mode function 

In the de Sitter background a = —l/(Hrj), the equa- 
tion of motion (EoM) of the mode function takes the form 

m, 



2 7 2 / 2 , r n 1 2 2 , u 7 4 4 

3 W 



(4.1) 

, (4-2) 



2 Quantization of the gauge invariant perturbation 1Z was per- 
formed in [l9j through the canonically normalized field v = 
z1Z,z 2 = a 2 h 2 (3o. For the simplicity of calculations, here we 
introduce the mode functions of 7Z, Tfc(??), which relates to Vk{rj) 
through the relations, 



rk(v) K 



(3.17) 



where c 2 s = bi and 62,63 are defined in [19]. Due to 
the time-dependence of the dispersion relation, the mode 
function Vk(rj) (or equivalently rk(t)) will not take a sim- 
ple plane wave form as in GR. Rather, its form will evolve 
with time. This complicates the calculations of the bis- 
pectrum considerably. 

One may worry that, like in the case of [16], relevant 
scales may have left the horizon at a time when the A: 6 
term dominated the dispersion relation. However, this 
seems not reasonable in the slow-roll inflation scenario. 
For the mode to leave the horizon at that time, two con- 
ditions have to be met, 

&34ifcV » <Z ~ h (UV regime) 
a" 2 

bse^k 6 T] 4 <C — = —5-, (Super- horizon region). (4.3) 
a r] A 

This would indicate that k 2 rj 2 <C 1 and at the same time 
bs ^> l/(e^ L /c 2 ?7 2 ) ^> 1. However, it was argued in [l9| 
that bs in general is of order one. Therefore, we consider 
this case as physically not realistic. 

Then, the evolution of the mode function during in- 
flation has to be taken into account. To deal with this 
problem, Brandenberger and Martin (BM) proposed a 
matching procedure [30j. Three regions were identified 
for the evolution history of the mode function (See Fig. 
[TJ More divisions are possible given a specific model.). 
Region I is the region in which the UV effects dominate, 
alias uj ~ k 3 in the present case. Then, the solution of 
the EoM in this region can be approximated with the 
Bessell functions of the first kind, 

vl(v) * Axy/\q\J v \z( V )} + A 2 y/h\J-» [<V)] ■ (4-4) 

In Region II, the dispersion relation restores its relativis- 
tic form uj ~ fc, so that the mode function can be safely 
approximated with plane wave solutions, 

Vk(rj) ~ Bi exp [—ic s krj] + B 2 exp [ic s krj] . (4.5) 

Note that unlike in the case of the Bunch-Davies vacuum, 
since the mode function underwent a UV stage, both pos- 
itive and negative frequencies appear in the mode func- 
tion in this region. Region III is the super horizon region 
when the mode freezes. The initial conditions in Region 
I are the ones that minimize the energy of the ground 
state of the field [30|, given, respectively, by Eqs. ([C.7|) 
and (|C.8|) in Appendix C of this paper, from which the 
two constants A\ and A2 are fixed. The undetermined 
coefficients of the solutions in Regions II and III are fixed 
by matching conditions across each boundary of these re- 
gions, by requiring that the mode function and its first 
order time derivative be continuous. The explicit expres- 
sions of these constants in terms of the initial conditions 
are given in Appendix C. 

On the other hand, the evaluation of the bispectrum 
all boils down to the following integration of the mode 
function, 

/** dt' 

rUt)rl 2 {t)rl 3 {t) J ^r fcl (t f )r k2 (t f )r ks (£')• (4.6) 



1 

H 


COph 

ti 

t 2 / 

















k 2 ki ki k 



FIG. 1: The evolution of cj p h = uJk/a vs k p h = k/a in three 
different regions, where Region I: £ G (ti, ti); Region II: £ G 
(£i,£2); an d Region III: £ G (£2,£o). 

A technical difficulty arises when the dispersion relation 
is of the form (j4.2ft as no exact solutions exist. Thus, 
the matching procedure presented in [30| seems a natural 
choice in approximating the solution of the mode func- 
tion to the EoM. Below we evaluate the mode integration 
(|3.16p and calculate the bispectrum and its shapes more 
explicitly using the BM matching procedure. 



B. The Bispectrum 

Dividing the integration into three regions, as men- 
tioned previously, one can see that only that over Regions 
I and II need to be considered, as Region III is the super 
horizon region, and the mode function gets frozen out. 
Then, we find that 

r* kl (tH 2 (t)r* k3 (t) I* ^r fel (t')r k2 (f )r fc , (f ) 
= U 1 (f e (ti, h); t) + U u (tf G (ti, t); t) , (4.7) 

where 

*W = r\\(t)r l l(t)T l l(t). (4.8) 

To carry out the above integrations, we first note 
that in Region L £ G (£i,£i), in contrast to the case 
considered in [16| where the modes all take the form 
exp[— i J dt'u/a 3 ] ~ exp[— ik 3 r]' 3 ], only those mode that 
are in the kernel of the integration take this asymptotic 
shape, while rj, 1 (£) takes the form of plane waves. There- 
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fore, in Region I we have, 



M 2 



U l (t'e(t h t 1 )-t) = — J2 B *v* E 



Iran 



x,y,z 
=1,2 



l,m,n 
= 1,2 



= K(k 1 )By(k 2 )B*(k 3 ) exp ( - ic s K xyz r) 
(J l lm n=Mki)A m (k2)A n (k 3 ) 
iH 2 



exp 



3M 2 



^imn fa? - 77?) 



(4.9) 



where the coefficients A and 5 are given in Appendix C, 
and 



K 



(-iyk( + (-i) m k* + (-i) n ki, 



5 = (-l)-^ + (-1)^ 2 + (-l)*fc 3 . (4.10) 



On the other hand, in Region II we find that 



x iVi z l.ra.n 
=1,2 =1,2 

\ rf -rjl 2 (77 — 771) 



II 

Iran 



i^trnn 



II _ 



crf^ n = Bi(k 1 )B rn (k 2 )B n (ks) exp ic s JCi mn (77 - 771) 

(4.11) 

Putting all the above together, we find that the bispec- 
trum can be cast in the form, 

(7e kl ft k2 7e k3 ) = 2/i 3 (2tt) 3 5 3 (kx + k 2 + k 3 ) 



x <6V, 



^1 ^2 ^3 



(fcj jfe| + k\k\ + cyclic) 



: 2 + ^3) 

xlmjf/ 1 + [/"}. 

(4.12) 



Using the expression of the mode function rk{t) given 
in Appendix C, we find that the non-linearity parameter 
/nl can be estimated as, 



7NL 



(Mill) 
(Ml) 2 



(27T) 3 



M* 



p2 
TZ 



X 



M 2 



(2ity 2tt 2 

H 2 ri\ (k + k + kf 



k 3 + k 3 + k 3 (k + k + kf 
xh 3 [QV 3 (k 2 k 2 k 2 ) + V 5 (3fc 6 - 6fc 6 )] 

° ( "^ 5) l~J (4.9 xlO-) 4 ' 



(4.13) 



where we have made the assumption k\ = k 2 = k 3 = k, 
/3q, c s ~ 1 and used the fiducial COBE normalization [27j 



of the power spectrum, from which we are able to repro- 
duce the spectrum presented in GR under the assump- 
tion H « M*. Writing M* = 10 _n M p i and assigning 
e ~ 10 -2 , when n = 2, we obtain 



10 



NL 



O (h 3 V 3j5 ) ~ 10 7 x O (h 3 V 3 , 5 ) ; (4.14) 



and when n = 3, we find that 

10- 7 



NL 



O (h 3 V 3 , 5 ) ~ 10- 1 O (h 3 ^ 3 , 5 ) . (4.15) 



Therefore, a large non-Gaussianity can be produced with 
a relative small V 3 and V5, given that the new scale M* 
isn't much lower than M p \. 



C. Shapes of the bispectrum 



Let us first note that, from (|4.12p . the /c-dependence 
of the spectrum, a.k.a the shape and running, depends 
not only on the action, but also on the form of the mode 
function. To study the /c-dependence of the bispectrum in 
more detail, let us first look at the expressions of U (£'; t) 
in Eqs. (|4.9j) and (j4.ll)) . The corresponding conformal 
time 77 is usually taken to be at very late of the inflation 
era, i.e., kn — >> 0. To be simple, here we take n = 0. 

With the coefficients A and B given in Appendix C, 
the bispectrum is plotted in Figures [2] - [5] for various 
choices of the parameters. Figures [2] and [4] are for the 
choice of parameters (|C.18p . whereas Figures [3] and [5] are 
for the choice of parameters (|C.19p . Within each figure, 
four sub-cases are plotted. 

From these figures, we can see that when S A and S B 
are taken to be zero, our result resembles the shape of a 2 
and a 3 terms in [10], despite the fact that our integra- 
tion (|4.7p is actually different from theirs. This is quite 
understandable from the following considerations: when 
we make S A = = 5 B , only the positive (or negative, 
depending on the choice of sign for the initial condition) 
frequency modes exist in both regions. This makes the 
product of the six mode functions in (|4.6p , which in gen- 
eral has 2 6 terms since each of them has two branches, 
collapses into one single term. Substituting this into Eqs. 
([HD and gTQ), we obtain 



U n (t' G (ti,t);t) oc 



1 



1 1 

(kik 2 k 3 ) 3 (ki + k 2 + k 3 ) 3 



(4.16) 



We see that the contribution from Region I has the same 
k- dependence as the a 2 and a 3 terms in (4.5) and (4.6) 
of jl_6|| . The contribution from Region II does not exist 
there, nor does it have the same ^-dependence as in the 
relativistic cases. 

The real difference comes in when we make either S A 
or S B nonzero. The impact of these on the shape of the 
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bispectrum is the enhancement of the folded shape [24[ 
(or sometimes called the flattened shape [28]), namely, 
that the bispectrum peaks at the limit k<i + ks — k\ — » 
and /cf + ^| — ^ 0. A non-zero 5 A for the positive 
frequency choice of the initial condition results in the 
appearance of "negative-frequency" modes in both r l ki (t) 
and (£), whereas a non-zero 5 B leads to a mixture only 
in r l ki (t). It's this mixture that makes the assumption 
which is usually kept in the standard choice of the Bunch- 
Davis (BD) vacuum 3 , that is, only positive frequency 
modes appear in the mode function, invalid. In fact, 
when we take rji to be infinite past like in the BD vacuum, 
S A will be zero automatically as can be seen from its 
definition in (|C.20|) . When we take S A = but 5 B ^ 0, 
this is similar to the case studied in [28| and our result 
is consistent with theirs. 

In addition, the choice of negative frequency A- es- 
sentially changes the sign of the bispectrum, and this is 
expected by comparing (|C.18|) with (|C.19|) 4 . 

We do not have an enhanced "squeezed" triangle signal 
as in the minimum scenario (though the overall magni- 
tude there is very small), nor do we have a large local 
form, which is a typical result of multi- field models. This 
is because the interaction terms that could generate the 
local form are all suppressed by factors of either c\jM v \ 
or e and are of sub-leading order [See discussions between 
(HTTP and (I3AQ1) ]. 

Now recall that in writing down (|3. 15j) . we have as- 
sumed that V s ' and V£ are zero. However, this is not 
physically necessary since Vs and V5 appears neither in 
the background equations, nor in the linear perturba- 
tions thus not constraint by the slow-roll conditions, and 
a strongly varying shape of Vs and V$ (not to be too 
strong to invalidate the perturbative expansion) would 
actually give both a higher non-linearity and new fea- 
tures in the bispectrum, such as the sinusoidal running 
and resonant running, as shown in [24j ]. 

Though the final integration was separated into two 
distinct periods, we would like to point out that this 
is a result of the matching procedure employed due to 
the lack of exact solutions for (|4.1|) in our model. One 
should, in principle, evaluate the integration as a whole 
and study the shapes of the bispectrum. A possible solu- 
tion is a development of the uniform approximation [34| 
with some numerical integration techniques involved. 



D. Projection onto factorizable templates 

To study the quantitative behavior of the shapes of 
the bispectrum, in general, factorable ansatz or template 



3 The initial conditions chosen in [l6j are the BD-like vacuum. 

4 Note that our self-interaction terms (|3.10|) have opposite sign 
w.r.t [l6|. This explains why our shape reproduces theirs only 
when we take the negative frequency. 




(a)S A = = 5 B 




1.0 

(b)S A = -0.1, 5 B = 




1.0 

(c)S A = -0.1, S B = 




1.0 



(d)S A = -0.1, 5 B = -0.1 

FIG. 2: Shape of the bispectrum (truncated). V3-term domi- 
nates, with the choice of the positive frequency. 
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FIG. 3: Shape of the bispectrum (truncated). V3-term domi- 
nates, with the choice of the negative frequency. 



FIG. 4: Shape of the bispectrum (truncated). Vs-term domi- 
nates, with the choice of the positive frequency. 
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(a)S A = = S E 




(b)6 A = -0.1, 5 B = 




(c)5 A = 0, <5 B = -0.1 




(d)(5 A = -0.1, 5b = -0.1 

FIG. 5: Shape of the bispectrum (truncated). Vs-term domi- 
nates, with the choice of the negative frequency. 



functions are utilized. Among them, three are of the most 
importance to us pij . 



Tbrth. (h,x,y) -18 



+18 



1 



1 1 



x 3 y 3 x 3 y 3 



-5—0 + 5 perms. , 

x z y 6 



48 
x 2 y 2 



T Fo i d .(k u x,y) = + \ + 



x 3 ?/ 3 x 3 ?/ 3 



18 

x 2 y 2 



x 2 y 3 



5 perms. , 



^Equi. (ki,x,y) 



,11 1 

-6 Mr + ^ 



+6 



1 



12 



x 2 y 3 



x^y* 

5 perms. I , (4.17) 



where x = k 2 jk\^y = k^/k\. Then the bispectrum ([4. 12)) 
can be cast in the form, 



(^ kl ^ k2 ^ k3 ) EE (2^) 3 £ 3 (ki + k 2 + k 3 ) 



xT H MT(*l,fc2,fc3), (4.18) 



where 



^HMT = CoTbrth. + PE^Fold. + CpTEqui. 



(4.19) 



with 



^HMT ' Tbrth. ^HMT " ^Fold. 

CO = — — , c F 



CE 



^Orth. • Tbrth. 
^HMT • ^Equi. 



^Fold. ' ^Fold. ' 



^Equi. ' ^Equi. 

where the inner product is defined as [l(| 

T 1 (k 1 MM)T 2 (k 1 MM) 



(4.20) 



rx-r 2 ee 

r-1 



oc x dx y dyTi(l,x,y)T 2 (l,x,y). 

JO Jl-x 

(4.21) 

When 6 A = = 5 B , the projection gives the following 
result 

c = M 2 ^1 + (0.0033T/ 3 - O.OOI8V5) ^, 

c F = M 2 ^1 + (0.0092T/3 - 0.0004V 5 ) h 3 , 
H 2 



c E = M [ 1 + ) (0.0063 V 3 - 0.0036 V 5 ) fc d , 



(4.22) 
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for the positive frequency, and 

c = M 2 ^1 + (-0.0004^3 + 0.0002^5) h\ 

c F = M 2 ^1 + (-O.OOOIV3 + O.OOOOF5) h 3 , 

c E = M 2 ^1 + (-O.OOO8V3 + 0.0005F 5 ) h 3 . 

(4.23) 

for the negative frequency. We see that the projection 
onto the equilateral template has the biggest magnitude 
for most of the cases. This is consistent with the obser- 
vations in [33j| that derivative coupling favors the corre- 
lation between modes with similar momenta. 

However, these analyses are only very rough approxi- 
mations. This is because these three bases do not form 
a complete basis set, as pointed out in [16], nor are they 
un-correlated, since we can compute correlations between 
them and get 

Corr (Equi., Orth.) - 0.204, 

Corr (Fold., Orth.) - -0.748, 

Corr (Fold., Equi.) - 0.489, (4.24) 

which indicate that there are strong correlations between 
the three, where 



Corr (SI, S2) 



(4.25) 



When 5 A and S B are not zero, the bispectrum, if taken 
naively, diverges at the folded limit /C2 + ks — k\ — >• 
or k 3 , + k 3 — k\ — >• 0, and the above inner product fails 
to converge. To understand the divergence, let us look 
at the integration (j4.9p in Region I more closely. The 
divergence occurs since, in addition to (j4.16p , we still 
have (2 6 — 2)-terms of the form, 

k3 + ls_ k s - (4-26) 

The denominator appears through the dv[ integration in 
(|4.7j) . When the denominator approaches zero in the 
folded limit, if S A and S B are not exactly zero, these terms 
blow up. The divergences must be regulated away if we 
want to have a quantitative estimate of the projection. 
A possible solution is to introduce a cut-off so that the 
integration yields, 

(O m (Ss) n ,3 ■ rs 1 ,3 I 5 (4-27) 



where 



k c (S A ,5 B ) , 



(4.28) 

in order to keep (|4.27p finite. As noted in [24], this 
regularization is highly model-dependent and requires 
some systematic study. Putting this divergence aside, 
we would expect that the projection onto the folded base 
would be enhanced significantly, in comparing with (|4.22[) 
and (|^23|) . 



V. CONCLUSIONS AND REMARKS 

We have studied in this paper non-Gaussianity of the 
primordial perturbations in a single scalar field with 
the slow-roll conditions in the framework of the Horava- 
Lifshitz gravity with the projectability condition and an 
arbitrary coupling constant A, refferred to as the HMT 
model. 

With some reasonable assumptions, we have found 
that the leading order terms in self-interaction of the 
HMT model are of order e 3 / 2 a n , where a n are dimen- 
sionless constants, defined as V n = a n /M^ v and V n are 
the coupling coefficients of sixth order derivative opera- 
tors of the scalar field, and have no contributions to the 
power spectra and indices, as shown explicitly in [l9| . 
Clearly, by properly choosing those coefficients, a large 
non-Gaussianity is possible. This is different from the 
standard result of minimum scenario in general relativ- 
ity, where the interaction terms are of order e 2 [zEUl- 

We have also investigated how the modified dispersion 
relation affects the evolution of the mode function and 
in turn the bispectrum, using the matching method pro- 
posed in 30] . By dividing the history of inflation into 
three regions, in which the dispersion relation takes dif- 
ferent asymptotic forms, the gauge- invariant scalar per- 
turbation 1Z has different asymptotic solutions. In partic- 
ular, we have found that the mode function is in general 
a superposition of oscillatory functions. This is different 
from the standard choice of the Bunch-Davis vacuum, 
where only one positive frequency branch of the plane 
wave is selected, and results in an enhancement of the 
folded shape in the bispectrum. 

Due to the existence of a UV region, the bispectrum is 
enhanced, and gives rise to a large nonlinearity parameter 
/nl [cf. Eq. (|4.13p ], as long as M*, above which the non- 
linear terms in the dispersion relation dominate, isn't 
much lower than the Planck scale [See (|4.14p and (|4.15p ]. 

We would also like to summarize the assumptions made 
along the analysis, as the invalidation of any of these 
assumptions will certainly change some of the conclu- 
sions. In particular, in obtaining the leading order term 
in the self-interaction, we made the assumption that 
c\ ^ M* <C M p i. This eliminates a great number of terms 
which are not of order Sx 3 in the third order expansion, 
and makes all terms that will give non-local effect (or 
local shape bispectrum) of sub- leading order. However, 
the phenomenological upper bound on M* for the pro- 
tectable version of the HL theory is not established yet 
[14l ]. Dropping this assumption will bring back a lot of 
terms and enhance the local shape signal in the bispec- 
trum. 

We've also kept the assumption made in [l9[ that ^ 
0. For Vi,^2 and V4, this assumption further eliminate 
some possible contributions to the leading order effect. 
For V3 and V5, as noted in the discussion at the end of 
Part.C in Section IV, if this assumption is dropped, new 
features (sharp and periodic shapes) will appear. The 
condition H/M* « 1 has been also assumed to obtain 
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an expression of the power spectrum similar to that of 
GR. However, if this is dropped out, the power spectrum 
will receive large corrections, as noted in [l9| . 

By dividing the time of interest into three regions, we 
have assumed that the period of dominance of the k A 
term in the dispersion relation is so short that its effect 
can be incorporated into the small parameter S B . This 
is a strong assumption, as noted in [35|, 36]. During this 
period of time, uj 2 may actually go below H 2 and then 
the mode function will be no longer oscillating, but grow 
with the scale factor a. This is an important problem 
and deserves further investigations. 
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Appendix A: Scalar perturbations in HMT Model 

For the reader's convenience, we present here the lin- 
earized perturbations and the corresponding constraints 
and field equations that were first obtained in [l9j . 
These will be used in obtaining the cubic action of self- 
interaction. Below we shall choose the Newtonian quasi- 
longitudinal gauge (j3.3p 5 . The variable ip defined as part 
of the perturbations to the 3-metric gij in that paper is 
related to ( used here via the relation, 

tl> ~ -(. (A.l) 

The constraint equations to first order are then given by, 



d 3 x\ d 2 ip - ^ (3A - l)% 3V>' + d 2 B 



= 4ttG 



J d 3 xffx'6 X ' + (a 2 V + §5 4 )^|, (A.2) 

(3A - 1)V' - (1 - X)d 2 B = SnGfx'Sx, (A.3) 
+ (l - A) (3*0' + d 2 B) 

= SttG [{c[x' + c x U - fx')S X + ci5 X '] , (A.4) 

ip = A-kGcxSx- (A.5) 



5 Note that the two popular gauges in GR: 

i) E = = 5x, B / 0, ^ and C + 0; 

ii) E = = C, B / 0, <j> / and 5 X / 

are not possible in our model because of (|A.5|) . 



The trace- and traceless-dynamical equations are, respec- 
tively, given by 



1 



i)" + m$ + \d 2 {& + 2UB) 



3(3A - 1) 
2 



^(fx'Sx' -a 2 V'5x), 



_8ttG 
3A~ 



where 



tp — (B' + 2HB) + 1 + ^d 2 )d 2 i> 
- i (Ai/> - SAj = 0, 



_ 852 + 3fi(3 _ 8g 7 - 3g s 

Ct\ = — 2 , «2 



(A.6) 



(A.7) 



M p 2 ,/2 ' M p V4 
The modified Klein-Gordon equation reads 

f{5 X " + 2H8 X ' - x! W + d 2 B] } + a 2 V"S X 

'l „ v* + vi^ v,\ a2 



(A.8) 



1 



d 2 



2A(c[ -c 2 )Sx + c 1 SA 



(A.9) 



Solving the above constraints and field equations yields 
B and other variables in terms of £, 



B 

SA 

a 



1 



(d~ 2 C) 



2/X' 



(1 - A)ci 



(a- 2 c) , (A.10) 



Vx 



12 



2V 



(1 — X)c\ c\ a 2 (l — A)ci 

2fx' __m 



(i-A)d 141 



Sx 



(ATrGd)-^, 



(A.11) 
(A.12) 



where c\ = (A - 1)/(1 - 3A), and (d~ 2 d 2 ) C = £. 
The quantity 



K 



v + —Sx, 

X 



( 1+ 47rGci^) C ' 



(A.13) 
(A.14) 



often referred to as the comoving curvature perturbation, 
is also gauge- invariant. It can be shown that in terms of 
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this quantity, the free action can be written as 

= i J dj]d z xa 2 h 2 \p 1l' 2 - (3 A TZ 2 - /3i(di1l) 2 
- h{d 2 U) 2 - fc{did 2 n) 2 



(A.15) 



where 



A) = f + ^Gc 2 /\c%\, 



/3i = l + 2V 1 + 2A(c' 1 -c 2 )-iTrGc 2 (l-A), 
= \(V 2 + Vl + 2TrGc 2 ai ), 

ft = (V 6 + 27rGc 2 a 2 ) , 



ft 2 



2 /i 2 



Q 



J_l 2t „, , 4TrGc lC '{x' 2 8ttG 



j/x ,2 (/- Cl ') 



-47rGcia 2 y' I 3 



f\c 2 



1 

l c V>l 



4» C e 1 + f)- =|. 



After introducing the normalized variable 



zTZ, z 2 = a 2 h 2 /3 , 



(A.16) 



(A.17) 



the classical equation of motion for the mode function v k 
takes the form, 



where 



v k + ("1 + m eff) ^fe = 0, 



Po 



-™eff 



Z Po' 



(A.18) 

(A.19) 
(A.20) 



Appendix B: Self interaction of the inflaton 

Under the Newtonian quasi- longitudinal gauge (|3.3p , 
with perturbations given by (j3.4p , we find that the action 
can be written in the form, 



^1(3) - Sglfs) + Sxlfs) 
+ S g\f3) + ^1(3) 



<7 I HL 



(B.l) 



where the "GR parts" are given by 6 
1 -3A 



C |GR 

^1(3) 



WttG 



J dr]dxa 2 i^Y^ 2 C 3 



9H( 2 (3C - d 2 B) - 2£ (d k Bd k () 



3C (C 7 ) - 2HC (d k Bd k () - 2(('d 2 B 

cm) 2 -Ac(^) 2 



1 - 


3A 


1-1 


-A 


1 - 

s 


3A 


1 - 


3A 



1 - 3A 
2 (d 2 B) (d k Bd k C) 



~\kc J d ^ dxa2 { 9Aa ^ 3 + C 292 C}, (B.2) 



C |GR 
^Xl(3) 



/ 



2 j 9c 3 r/(A) 



drjdxa < 



X ,2 -a 2 V 



9^ 
2 
3C 



f(\)x'5 X ' - a 2 V'5 X 

2 



/(A) (5 X 'Y - a 2 V"5 X 2 
3Cf (A)x' (d k B) (d k S X ) 
f(X)5 X ' (d k B) (d k 5 X ) 



(B.3) 



And the "HL gravitational part" is given by, 

^ J d V dxa 2 f -a 2 J? v \™ 

(16 5 2 + 5 ff3 )(9 2 C) 5 



9 '(3) " 16ttG 



(30 



a 2 M 2 



■93(dij0 
2_\ 2 



5# 8 (d k d 2 () +g 8 (d ijk CY 



(B.4) 



Note that though this is labeled as the "GR part", it cannot 
reproduce the exact expression of GR, as a result of the difference 
in symmetry. 
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where 



(_ 6 4 52 -2o 53 )C(a 2 c) 2 



+ (16 52 + 6 53 ) (d 2 C) (dkC) 2 

-2 ff3 (^C3 i C^C)-4 5 3C(%C) 2 
+ 16 58 (diCdjd^d^c) 

- (32 57 + 20 58 ) (5 2 Ca fe C9 fe a 2 C) 

-30 58 c (a fe a 2 c) 2 -6 58 c(a ljfe c) 2 

- (64 54 + 20 55 + 6 56 + 32 57 ) (9 2 C) 3 
+ (1657 - 4g 5 - 3 56 ) (9 2 C) (d tj Cf 

-56 (a^kC^C) 

- 96 57 ca 2 ca 4 c + 8g 7 d 4 c (dkO 2 ' 



(B.5) 



The "gauge part" is given by, 



5a |(3) = I dr 1 dxa 2 {2A{ 2 (d 2 <;)+2A<;(d k 2 



4((d 2 ()SA + 2(d k 2 5A\. (B.6) 



Finally, the "HL matter part" is given by 



q |HL 



/ 



drjdxa 



2 j9C 2 



v 4 



4 (a 2 ^) 



+ (3a 2 C) JS? X |$ +a 2 



(B.7) 



|HL 



2a" 



.^2 



_ 2 _± (^x) [7>i| (2) 



YL 

a 4 



(<Jx) (9 2 ^) 2 
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^2 1 (3) 



-5 (^) 8 

-Vk'(*X) [^2 1(2) 

~{^S X ) (d k S X ) 2 

-§(*x) (d 2 s x ) (dH x ) 
v 6 



^4 Vxf (3%) 



I c'/ 



a 2 2 



(0 2 «J X ) (a 4 P 2 | (2 ))+(^x) (a 2 Pi| (2) ) 
Pi | (3) ] + Ac' 1 (8 X ) [Vi\(2) 

(Sxf (d 2 x) 

2c 2 { + c' 2 6x] (d k S X ) 2 



a 



c[S X (d 2 5 X )+c 1 a 2 r 1 \ {2) 
+ c 2 (d k S X ) 2 ], (B.9) 
where 

a 2 Pi| (2) = (dk() (d k 5 X ) - 2C (d 2 5 X ) , 

a 4 P 2 | (2) = [(d k 5x)(d k d 2 ()+2(d k O(d k d 2 Sx)' 

+2(^-0 (^'5x)-4C (9%), 
a 2 Pi|( 3) = 3C 2 5 2 5 X , 
« 4 P 2 |( 3) = $C 2 d 4 5 X - 8C (%C) (0*%) 

-4C [ (d k 5 x ) (d k d 2 () + 2 (d k C) (d k d 2 5 x ) 

+ (dX) [ {d j C) (dijSx) + i&Sx) (%C) 

(B.10) 

One can further substitute 5y, B and 6 A in terms of 
C into Eqs. (fA7TO|) - ([A7T2]) . See [37| for a very good and 
detailed review of the subject. 



with 



|HL 

^xl(2) 



"V 4 



^2 | (2) 



^ + ^i)a- 2 (^ X ) 2 



a 

cxS A 



Pi I (2) 



^ C2 ^ JC \2 



(fl 2 *x) , 



(B.8) 



Appendix C: Determination of the integration 
constants 

We find that in the original work [30|, one of B\ and 
B2 must be zero in the case of the Corley-Jacobson dis- 
persion with bm > 0. In fact, considering their equations 
(134) and (135), we find that the terms in the parentheses 



B 



1 oc 



(cos 2/1 



1 sin yi 



e^cos x x 



(1 



^2 



oc (cos yi 

= e*' y± — e n 



1 sin 7/1 



e Ti7ri/ cos x x 



-ze^^sin xx) 



■ ze Ti7ri/ sin xi) 



(C.l) 
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where we have used the definition of x and y in equa- 
tion (128), x = y + 7rv. One immediately sees that either 
B\ ^ 0, B2 =0 (upper signs, no excitation and only "pos- 
itive frequencies" exist), or B\ = 0, B2 7^ (lower signs, 
complete excitation and only "negative frequencies" ex- 
ist). Similar problems occur for their A± and A2 too 
(take the large argument asymptotes of the Bessell func- 
tions). We believe this is an artifact of the approxima- 
tion procedure, and that in general, both the "positive-" 
and "negative-frequency" modes should appear and the 
magnitude of the "negative-frequency" modes should be 
limited by considerations from the back-reaction problem 

mum. 

Below we re-derive the coefficients following the BM 
guidelines of matching, but within our concrete model. 
From (|A.15[) . the mode function rk(t) satisfies the EoM 



3H 



rn 



a 2 



r*(t) =0, 



(C.2) 



where m 2 = /?4//3o> w\ is the same as defined in Eq. 
(|A.19p and a dot represents differential w.r.t time t. In 
Region I where the nonlinear effects in the dispersion 
dominates, the solution takes the asymptotic form (below 
we use t and conformal time rj interchangeably, noting 
that adrj = dt.), 



r\(t) 



Vi 



A! 

A 2 

M" 4 



■ H 2 k 3 ( o 
.H 2 k 3 ( o ,\ 



A2 exp 



A2exP L ~3M 2 
fta 4 _ 2 1^ + 2- 
~ft~ ~ 




2 u3 



.H 2 k 
' 3M 2 



f + ^Gc 2 /\cl\ 



(C.3) 
(C.4) 



It can be shown that M here is related to M* defined in 
HH through Mt = b 3 M 4 . Since 63 ~ 0(1), the condition 
H 2 /M 2 < 1 implies H 2 /M 2 <C 1. 

In Region II, the mode function is a superposition of 
plane waves, 



Bi 
B 2 

2 



Bi exp [-ic 3 k (rj - rji)] + B 2 exp [ic s k (rj - 771)] 
B\ exp [—ic s kn] + B 2 exp [+ic s krj\ , 
B\ exp [+ic s krji] , 

B 2 exp [-ic s krji] , (C.5) 
ft _ 1 + 2Vj + 2A(c' 1 - c 2 ) - 4irGc 2 (1 - A) 
ft" / + 47rGc 2 /|c 2 | 

(C.6) 



The initial conditions are chosen such that 

Vk(Vi) 1 



rk(U) 



Km) z^y/MmY 



(C.7) 



whereas its initial time-derivative takes value such that 
the energy density is minimized at the initial time (not 
necessarily infinite past rji — » — oo) 



= ±i- 



(C.8) 



z(Vi) V 2 

At the time of matching between region I and II, ti, 

r l k (t 1 )=r%(t 1 ), 4{t 1 )=4\t 1 ). (C.9) 

These four conditions ()C.7|) ^ (jC.9|) fix the four undeter- 
mined constants 



exp 



• H k 3 
1 3M 2 '/i 



:(77i) (77*) V 



(> 



2i 



exp 



3M. 







2z 



M 2 



(CIO) 



for r' k (U) = 



(positive frequency choice), 



A 2 



;H z k s ^3 



exp ^cI-zmT'H / 2Z \ 

- 1 — * J + , 

7i) V^cj (rji) v r]i^(m) J 



exp 



3M2 '/i 



for r^) 
And 

£1 



^7) 



( 



2z 



).(C.ll) 



-i—^y (negative frequency choice). 



a +ic s kr)i 



k 2 H 2 rjl\ ~ ( H 2 k 3 , 
1 + ttt^- I ^1 exp ( -t^j^Vi 



c s M 2 



3M2 



k 2 H 2 ril\ ~ ( .H 2 k 3 , 
1 rrrf 1 I A 2 exp +z 



c.M, 2 



3M2 



»7i 



B 2 



-ic s krji 



k 2 H 2 V 2 \ ~ ( H 2 k 3 , 
1 rr^- I ^1 exp ( —t-^Y7 2 ~ r li 



CsM 2 



3M2 



f^)i 2 exp +l - 



c s M 2 



2^1 



(C.12) 



The time of matching 771 (or equivalently, ti) is a critical 
quantity, here we choose it to be [32| 



^M-c 2 fc 2 = c 2 fc 2 , 



(C.13) 
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or from ([A7L9]) 



or 



(C.14) 



Now recall that by dividing the evolution history into 
three regions, we implicitly assumed that the /?2 term 
was never dominant during the history of inflation, as a 
result, the above condition can be approximated with 



which leads to 



2 \V + 2irGcla 2 \ k^H 4 ^ = ft. 



(C.15) 



(C.16) 



Hence 



(kH m f = c s M, 2 (14 
711 * - — ~H 



2<5 B ), 
(1 + *b), 



(C.17) 



where we have used the definition of c s and M*, and 
introduced a quantity S B to denote any deviations from 
(|C.15p . for example, when we consider the k A term's min- 
imal effect on 771. (In [32], it was found that a difference 
in this matching time could result in an unnecessary os- 
cillatory component for the power spectrum, and we shall 
see below that this indeed happens through an extra os- 
cillatory phase for the coefficients B.) 

Looking at the equation (|C.12|h this choice of 771 indi- 
cates that the influence of A2 (Ai) on B\ (B2) is minimal. 
Also note that rjiUJ (rji) ~ (H 2 /M 2 )k 3 r]f is in general a 
very large quantity, this allows us to make further ap- 
proximations on the expressions of the coefficients A and 
B and arrive at 



A 2+ 
B 2 + 



—M* 

hv%V2k? 
—M* 



e +i4>A (1 + 15a) > 

e-* A (i5 A ) , 

e^ B1 (1 + i5 A + S B ) , 

(i<5Ae- #B1 +5 B e i,pB2 ) , 



(C.18) 



A 2 - 
Bi- 
B 2 - 



where 



—M* 
M, 
—M* 



(iS A ) . 



e- i<tfA (l+iS A ), 
(iS A e i(f)B1 + 5 B e~ i(f)B2 ) , 
e -^ B1 (l+iS A + S B ), (C.19) 



2M 2 



rji(j(rji) k 3 r] 3 H 2 ' 
H 2 k 3 , 



(C.20) 



3M + 2 

c s kr]i 
3 

Cskrj! 



Vi 



(1 + 25 B )(^) +2(1-J fl ) 



(1 + 2S B )[^) -2(2 + <J B ) 



(C.21) 
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